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Today’s talk: CFD applications to launch environment focusing on 

unsteady flow in flame trench and prediction of ignition 
overpressure 
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Outline 



• Introduction / Background 
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- STS-1 and STS-1 24 Analysis 

- Ares- IX 
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Introduction: Major Sources of Unsteady Loads on Launch Systems 


1. Ignition Over Pressure (IOP) 

Primary source of transient load on flame trench, 
launch structure and vehicle at launch 

• High-accuracy simulation of start-up process was first 
performed at Ames in conjunction with Ares-IX 

2. Plume aero-acoustics 

Can impose significant vibration load on vehicle and 
payload during ascent 

• Experimental correlation based on a single rocket is 
being used: NASA SP-8072 (1971) 

3. Vibration due to turbopump 

Can generate damaging vibration on engine and the 
whole launch vehicle during ignition, lift-off and ascent 

• Prediction tool not available to quantify cavitation 

4. Combustion-related vibrations, such as thrust 
oscillation and combustion instability 



IOP 



Sketch of overall sound pressure 
contour from exhaust plume 



SSME Low Pressure Fuel Pump: instantaneous 
pressure map & measured pressure fluctuation 
upstream of the pump 
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Introduction: Unsteady Flow in Flame Trench 


IOP Wave Physics 

- During ignition, the exhaust 
plume injects mass into the 
confined volume of the flame 
trench under the Mobile Launch 
Platform (MLP). 

- This additional mass displaces 
the air within the trench causing a 
piston-like action. 

- Compression waves then travel 
up and down the trench 
generating a series of strong 
pressure waves. 

- The pressure waves travel back 
through the MLP exhaust holes 
towards the launch vehicle, 
possibly damaging the vehicle 
and surrounding structure 




Shuttle and flame trench: main definition 



Side view 
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Introduction: Unsteady Flow in Flame Trench 



In addition to IOP wave phenomenon, other complex physical 
processes also occur during ignition and liftoff. 

- Transient build up of the chamber stagnation conditions (very 
fast time scales). 

- Multispecies interaction of exhaust gases with the ambient air. 

- Fuel rich exhaust gas afterburning (chemical reactions). 

- After burning effects of solid aluminum particles in SRB 
plumes. 

- Multiphase interaction of the exhaust gases with the IOP wave 
suppression system (for example, water jet, water bag etc). 
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Introduction: Historical Background 



STS-1 & 2 Flight Data Comparison 

- STS-1 data shows very high ignition 
pulse 

- STS-2 shows the impact of water 
suppression system 



S j 

CFD simulation of IOP with 
Shuttle, Kiris et al. 

STS-1 April 12 , 1981 
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STS-1&2 IOP Flight Data 
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Introduction: Historical Background 



STS-124 Trench Wall Damage 

During the launch of space shuttle Discovery 
on May 31 , 2008, flame trench wall was 
damaged. 

Debris shown in the photo is the residue from 
this damage. 

Repairs are done before space shuttle 
Atlantis’ STS-125 mission to NASA’s Hubble 
Telescope. 



Debris scattered near 
Launch Pad 39A at NASA's 
Kennedy Space Center after 
May 31 , 2008 launch of 
space shuttle Discovery 




Damaged wall 


Debris scattered outside the perimeter of 
Launch Pad 39A 
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Prediction of Ignition Over Pressure (IOP) Waves 



• Objectives of CFD Simulation 

- Characterize/ quantify IOP wave for new and existing launch vehicles 
(as a part of the entire ascent simulation). 

> IOP trend analysis can contribute to launch pad design. 

• Approach 

- Evaluate current capabilities and develop high-fidelity methods/tools 
Algorithm and solution procedure 

Space time resolution requirements 

- Idealized test case 

2-D impinging jet: Experiment by JAXA 

- STS-1 case 

- Single-phase applications (Ares-IX) 

- Multi-phase modeling issue (STS-4) 


(Currently collaborating with JAXA to develop numerical methods, 
enhance supercomputing performance, and develop advanced multi- 
phase flow modeling) 
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CFD Procedures 



• State-of-the-art CFD technology 

- Grid generators, solvers (and associated algorithms) are generally 
available (either in-house or commercial codes). 

- Engineering-level physical models are also available for steady state 
solutions. 

- Prediction capability is limited to small regions, primarily in steady- 
state (generally good for interpolation). 

- Supercomputers at Petaflops level are becoming more available, 
requiring to look into parallel efficiency, data management... 

Jaguar at ORNL (224,000 cores) 

Pleiades at NASA Ames (81,000 cores, theoretical peak of 973. 3TF) 

...and others 

• Current Procedure 

- Rigorous error estimate/control methods, validation procedures are 
not available - “Best Practices” are current approach. 

> Next, will demonstrate how space and time resolution requirements 
are determined for IOP simulation. 
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Computational Model and Flow Solver 



• Computational Model and Grid: 

Overset grid: 

- e.g. STS-1 Shuttle configuration 

• 2 SRB’s and external tank 

• 129 overset grids 

• 92 Million grid points 

• Flow Solver : OVERFLOW-2 

- Compressible Navier-Stokes flow solver 

- Automatic domain connectivity including 
grid splitting for parallel computations 

- Spalart-Allmaras Turbulence Model 

- Diagonalized Beam-Warming scalar 
pentadiagonal scheme for LHS terms 

- Second-order in time with sub-iteration 
procedure (20 subiterations) 

- Physical time step used : 8.0e-5 seconds 

• Assumptions 

- Single species calculations 

- Single Phase (no water effects) 
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Governing Equations 



Reynolds Averaged Navier-Stokes Equations 
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Governing Equations 
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Time Step Sensitivity Analysis 



Dual Time Step Formulation in OVERFLOW 2 

dQ + d{F-F v ) | d(G-G v ) [ <?(£-£„) Q 

t?T dt dfj d £ 


Generalized Coordinate Transformation 
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Time Step Sensitivity Analysis 


Discretizations can be done, e.g. by fully-discrete formulation 
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Time Step Sensitivity Analysis 


Test Problem-2D Flame Trench 

Time-Accurate Pressure Point Locations 



Initial Conditions: Atmospheric Pressure and Density, Zero Velocity. 


Boundary Conditions: Time-accurate plenum conditions at the nozzle. 
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Time Step Sensitivity Analysis 
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Time Step Sensitivity Analysis 
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Time Step Sensitivity Analysis 
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Time Step Sensitivity Analysis 


Analysis is performed using the dual time stepping framework to analyze the 
sensitivity of unsteady IOP wave propagation with respect to both time step and 
convergence. 


• The largest possible pseudo time CFL numbers should be used for efficiency 

• Larger pseudo time CFL numbers were examined, but lead to sub-iteration 
instabilities. 

• Large physical time steps with small pseudo time CFL numbers and a small 
number of sub-iterations lead to nonphysical results (without numerical 
instability). 

• For a fixed time step and pseudo time CFL, the sub-iterations should be 
increased until a sub-iteration invariant solution is obtained. 


1 ,0e - 05 ^ At ^ 5 .0e - 05 for nsubs = 10 and 20 
5 .0e - 05 ^ At < 1 .0e - 04 for nsubs > 50 
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Idealized Test Case for Validation 


Test Case Introduced by N ASA- JAXA collaboration: 

2D Impinging Jet 

• Focus on basic methods for capturing IOP and launch 
acoustics 

• Conical Mach 2 nozzle 

• Single-species (y = 1 .4) 

• Po/Patm = 7.825 

• Over inclined flat plate 



M d =2.0 


Source: JAXA 
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Simulating Ignition Overpressure Waves 


Numerical experiments 

•Jet Impingement Physics 
•Sensitivity Analysis 
•Sub-Iteration 
•Physical Time-Step 
•Space-Time 




Parameter Matrix 

- Seven sub-iteration counts: 

• NSUB = 5, 10,20,50, 
100, 200,400 

- Six physical time levels: 

• CFL= 1, 2, 5, 10, 50, 
100 

- Three grid resolution levels: 

• Coarse Ax = 0.005 m, 

1 .4e-05 <= At <= 1 .4e-04 

• Medium Ax/2, At/2 

• Fine Ax/4, At/4 


. CFLAx 
A t = 

Cref 
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Jet Impingement Time Sequence 
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Gauge Pressure (psig) 


Sub-Iteration Sensitivity Test 



Location 6 and 10: 
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Time-step Sensitivity Analysis 



Location 6 and 10: 
o CFL <= 10: Appear Converged 
o NSUB/CFL >= 40: Sufficient 

Overset 

Location: 6 Grid: Medium NSUB: 0400 



Location: 10 Grid: Medium NSUB: 0400 
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Space-Time Sensitivity Analysis 


Location 6 



Coarse Medium Fine 

• Low-frequency * Increasing frequency • High frequency 

develops after IOP 


29 


Idealized Test Case 



• Numerical experiments have been performed for an idealized test problem 
Comparison with experiments will be done as data become available 

• Preliminary conclusion from the numerical tests: 

- For a fixed ratio NSUB/CFL, the sub-iteration convergence is approximately 
independent. 

- If time-integrated functionals (e.g. RMS) are of interest, then coarser space- 
time resolutions may be acceptable. 

-When high-frequency wave content is important (e.g. Min/Max functionals), 
then viscous effects should be included. 
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Applications to Shuttle and Ares-IX 



Assumptions for CFD simulations 

- Single species flow with corrected nozzle boundary 
conditions such the correct thrust, temperature, and 
Mach number are retained at the nozzle exit. 

- Influence of multispecies gas effects are negligible. 

- Influence of solid aluminum particles in the exhaust gas 
and afterburning effects are negligible. 

- IOP water suppression system is neglected, thus no 
multiphase flow effects are included in the present results. 

- All relevant geometric structures are included in the model. 

- Full unsteady RANS simulation is carried out. 



Flow patterns showing fountain effects 
in impinging plumes 


Required information for CFD simulation 

- Time accurate conservative quantities for a perfect gas at the 
nozzle plenum. 

Quantities prescribed at the plenum must generate the correct 
thrust, temperature, and Mach number at the nozzle exit. 


• Given information 

- Pressure and mass flow rate at the nozzle plenum. 

- Mixture exhaust gas properties. 

- Approximate temperature of exhaust gas at the plenum. 
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Flame Trench Flow Analysis for STS-1 
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STS-1 IOP Flight Data vs. Computations 
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Flame Trench Flow Analysis for Ares-IX 




Pressure Contours 
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Flame Trench Flow Analysis for Ares-IX 




Pressure Contours 


35 



Flame Trench Flow Analysis for Ares-IX 




Pressure Contours 
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Pressure 


STS-1 andAres-IX IOP computations 



IOP Inner Side of SRB above nozzle skirt 



Time (seconds) 


37 


CFD Support for STS-125 Launch Environment 


After STS-124 wall damage on May 31 , 2008, trench wall had to be repaired 
before STS-125 mission to Hubble Telescope. 


• To quantify flow in flame trench (pressure, 
temperature etc.) 

- Computations were carried out up to 1 .1 5 
seconds by using 504 CPU’s (completed 
in less than 4 days) 

- CFD results were compared against 
STS-4 flight data. 

• Single-phase computations produced 
conservative results compared to STS-4 
data. 

• Time-dependent computational data were 
provided to Ground Operations to determine 
load environment. 
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Computed Results (Overflow) vs. STS-4 Data 
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Computed Results (Overflow) vs. STS-4 Data @ f 
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Computed Flame Trench Wall Pressure 



Time: 0.00000 



Computed Flame Trench Flowfield 
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Summary and Discussions 


Computational processes/issues for supporting mission tasks are discussed using an 
example from launch environment simulation 

• Entire CFD process has been discussed using an existing code 

• STS-124 conditions were revisited to support wall repair effort for STS-125 flight 

• When water bags were not included, computed results indicate that IOP waves with the peak 
values have been reflected from SRB’s own exhaust hole. 

• ARES-1X simulations show that there is a shock wave going through the unused exhaust hole, 
however, plays a secondary role 

• All three ARES-1X cases and STS-1 simulations showed very similar IOP magnitudes and patters 
on the vehicle. With the addition of water bags and water injection, it will further diminish the IOP 
effects. 

For more complete flame trench simulation, need to include 

• Water suppression system 

• More complete nozzle (and plum) condition 

• Multi-species (variable y) effects 


For more predictive launch environment analysis (supporting space exploration system 
development and operation in general), need 

• Predictive unsteady flow simulation capability (algorithm + physical model) 

(This is a current CFD pacing item) 
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